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Abstract 

The dynamical density fluctuations around the QCD critical point (CP) are analyzed using rel- 
ativistic dissipative fluid dynamics, and we show that the sound mode around the QCD CP is 
strongly attenuated whereas the thermal fluctuation stands out there. We speculate that if pos- 
sible suppression or disappearance of a Mach cone, which seems to be created by the partonic 
jets at RHIC, is observed as the incident energy of the heavy-ion collisions is decreased, it can 
be a signal of the existence of the QCD CP. We have presented the Israel-Stewart type fluid 
dynamic equations that are derived rigorously on the basis of the (dynamical) renormalization 
group method in the second part of the talk, which we omit here because of a lack of space. 



1. Introduction 

A unique feature of the QCD phase diagram is the existence of a critical point. At the QCD 
CP, the first order phase transition terminates and turns to a second order phase transition. Around 
a critical point of a second order transition, we can expect large fluctuations of various quantities, 
and more importantly there should exist a soft mode associated to the CP. The QCD CP belongs 
to the same universality class as the liquid-gas phase transition point, and, hence, the density 
fluctuating mode in the space-like region is a softening mode at the CP: The would-be soft mode 
of the chiral transition, the cr mode, is coupled to the density fluctuation! 1 ] an d becomes a slaving 
mode of the density variablejj]; see yl] for another argument on the fate of the cr mode around 
the CP. 

The density fluctuation depends on the transport as well as thermodynamic quantities that 
show an anomalous behavior around the critical point. In particular, we should note that the 
density-temperature coupling which was not explicitly taken into account can be important. In 
fact, the dynamical density fluctuations are analyzed in the non-relativistic case with use of the 
Navier-Stokes equation, which shows that the Rayleigh peak due to the thermal fluctuation would 
overwhelm the Brillouin peak due to the sound modes||4j]. 

We apply for the first time relativistic fluid dynamic equations to analyze the spectral prop- 
erties of density fluctuations, and examine possible critical phenomena. We shall show that even 
the so called first-order relativistic fluid dynamic equations have generically no problem to de- 
scribe fluid dynamical phenomena with long wave lengths contrary to naive expectation. In this 
report(5|], we shall show that the genuine and remaining soft mode at the QCD CP is not a sound 
mode but the diffusive thermal mode that is coupled to the sound mode, and that the possible 
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divergent behavior of the viscosities might not be observed through the density fluctuations be- 
cause the sound modes are attenuated around the CP and would eventually almost die out at the 
CP. 



2. Relativistic fluid dynamic equations for a viscous system 

The fluid dynamic equations are the balance equations for energy-momentum and particle 
number, <9 Jli 7 vv = 0,8^** = 0, where is the energy-momentum tensor and the particle 
current, respectively. They are expressed as T^ v — (e + P)u^u v - Pg^ + t' jv and ffl = nu^ + yP, 
where e is the energy density, P the pressure, uf the flow velocity, and n the particle density, the 
dissipative part of the energy-momentum tensor and the particle current are denoted by r /iv and 
v^, respectively. 

The so called first order equations such as Landau|0 and EckartJ3l equations are parabolic 
and formally violates the causality, and are hence called acausal. The causality problem is cir- 
cumvented in the Israel-Stewart equation|8], which is a second-order equation with relaxation 
times incorporated. One should, however, note that the problem of the causality is only encoun- 
tered when one tries to describe phenomena with small wave lengths beyond the valid region 
of the fluid dynamics: The phenomena which the fluid dynamics should describe are slowly 
varying ones with the wave lengths much larger than the mean free path. Indeed, the results for 
fluid dynamical modes with long wave lengths are qualitatively the same irrespective whether the 
second-order or first-order equations are used or not[5]. As for the instability seen in the Eckart 
equation|9], a new first-order equation in the particle frame constructed by Tsumura, Kunihiro 
and Ohnishi (TKO) jfioll has no such a pathological behavior. We employ Landau[6], Eckart|01, 
Israel-Stewart(I-S)|[lrand TKO equation. 



3. Spectral function of the dynamical density fluctuation 

By linearizing the fluid dynamic equation around the equilibrium, we can obtain the spec- 
tral function of the density fluctuation. The calculational procedure is an extension of the non- 
relativistic case described in the text book j4fl. 

The spectral function derived from the Landau equation is found to be 

S m (k,(o) = ((6n(k,t = 0)) 2 )[(l--) 2FrA " 



i r B fc 2 | r B k 2 
+ j (w - c s k) 2 + r 2 B k 4 + (to + c s k) 2 + rlk 4 ' " ' 

Here, the first factor represent the static spectral function, which would show a divergent behavior 
in the forward angle (k = 0) at the CP; this is known as the critical opalescence. The first term in 
the square bracket represents the thermal mode called Rayleigh mode, whereas the second and 
the third the sound mode or Brillouin mode. 

The Eckart equation in the particle frame does not give a sensible result for the dynamical 
density fluctuation, in accord with its pathological property[9). It is noteworthy that newly pro- 
posed equation, the TKO equation lUOIl . in the particle frame gives a sensible result even thou it is 
a first-order equation. We have also applied the Israel-Stewart equation!!! in the particle frame 
to obtain the spectral function for the dynamical density fluctuation. The result is the same as 
that of Landau equation; this tells us that the modified part to circumvent the causality problem 
does not affect the dynamics in the proper fluid dynamic regime. 
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4. Critical behavior of the dynamical density fluctuations 



We examine the critical behavior of the spectral function of the density fluctuations around 
the QCD CP. We introduce the static critical exponents y and a which are defined as follows 
c n - cot~ a , Kj = Kot~y, where t = \(T - T c )/T c \ is a reduced temperature, cq and Kq are constants 
and Kj = (1 / «o)(<9n / dP)j is the isothermal compressibility. We also denote the exponent of the 
thermal conductivity by a K , i.e., k = Kot~ a ', where kq is a constant. It is known that a K ~ 0.6 
around the liquid-gas phase transition point. 




-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 

o)[1/fm] ra[1/fm] 

Figure 1: The spectral function at t = (T — T c )/Tc = 0.5 (left panel) and at t = 0. 1 (right panel) for k = 0. l[l/fm]. The 
solid line represents the Landau/Israel-Stewart case, while the dashed line the TKO case. The strength of the Brillouin 
peaks due to the machanical sound mode becomes smaller as T approaches T c due to the singularity of the ratio of 
specific heats; the Brillouin peaks eventually die out as seen from the right panel, where the difference between the 
Landau and TKO cases is not seen anymore. Note that the scale of the vertical line in the right panle is much bigger than 
that of the left panel. 



Unfortunately or fortunately, these singular behaviors of the width of the Brillouin peaks 
around the QCD CP may not be observed. The strengths of the Rayleigh and the Brillouin peaks 
are given in terms of y as seen from eq.©, the ratio of the specific heats, which behaves like 
y = c p /c„ ~ r^ +a — > oo, in the critical region. Then the strength of the Brillouin peaks is 
attenuated and only the Rayleigh peak stands out in the critical region, as shown in Fig.Q] 

Let £ = gof~ v be the correlation length which diverges as the critical point is approached. 
If we write the wave length of the sound mode by A s , the fluid dynamic regime is expressed as 
£ << A s , with which condition the sound mode can develop. However, in the vicinity of the 
critical point, the correlation length £ becomes very large and eventually becomes infinity, so the 
above inequality can not be satisfied, and the sound mode can not be developed in the vicinity of 
the critical point. 

From this argument, we can speculate about the fate of the possible Mach cone formation 
ifTlll by the particle passing through the medium with a speed larger than the sound velocity c s . 
Such a Mach-cone like particle correlations are observed in the RHIC experiment llla . Then the 
disappearance or suppression of the Mach cone according to the lowering of the incident energy 
by RHIC would be a signal of the existence of the QCD critical point provided that the incident 
energy is large enough to make parton jetsQ5|] . 
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5. Concluding remarks 

In this report, the density fluctuations is analyzed using the relativistic fluid dynamic equations!]!] . 
We have suggested that a suppression or disappearance of the Mach cone formation with lower- 
ing the incident energy at RHIC can be a signal of the detection of the QCD CP. Although we 
have presented the Israel-Stewart type fluid dynamic equations that are derived rigorously on the 
basis of the (dynamical) renormalization group method, we omit them here because of a lack of 
space. For the details, we refer to the submitted paper B13ll . where it is shown that the transport co- 
efficients have no frame dependence while the relaxation times are generically frame-dependent 
in the derived equations. 
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